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Abstract. Let G be a reductive algebraic group of Q-rank one associated to a 
self-adjoint homogeneous cone denned over Q, and let T C G be a torsion-free 
arithmetic subgroup. Let d be the cohomological dimension of T. We present 
an algorithm to compute the action of the Hecke operators on H d (T: i 'L). This 
generalizes the classical modular symbol algorithm, when T C SX 2 (Z), to a setting 
including Bianchi groups and Hilbert modular groups. In addition, we generalize 
some results of Voronoi for real positive-definite quadratic forms to self- adjoint 
homogeneous cones of arbitrary Q-rank. 



1. Introduction 

1.1. Let $) 2 = SL 2 (R)/SO(2) be the upper-half plane, and let T(N) C SL 2 (Z) be 
the principal congruence subgroup of level N > 2, acting on Sj 2 from the left by 
linear fractional transformations. Then the cohomology group H 1 (V (N)\$)2] C) is 
closely related to the space of all weight-two modular forms of level N. The modular 
symbols provide a concrete approach to the group H 1 (T(N)\9)2', C) (§2.1]) that has 
allowed the testing of many conjectures in number theory and has led to explicit 
formulas for L-functions and their derivatives H [PH Hl8| ■ Important to applications 
is the modular symbol algorithm developed by Manin [18|]. An algebra of Hecke 



operators acts on H 1 (T(N)\Sj2', C), and using the algorithm one may compute their 
eigenvalues. Essentially this algorithm is the euclidean algorithm applied to pairs of 
integers (j p.2f) . 

Now consider the case where T is a torsion-free Bianchi subgroup, that is, T is of 
finite index in SL2(0k), where &k is the ring of integers in an imaginary quadratic 
extension K/Q. The group T acts on hyperbolic three-space $js = S X 2 (C) / 'SU "(2) , 
and we consider the cohomology group if 2 (T\f)3; C). As before, an algebra of Hecke 
operators acts on the cohomology, and one is interested in this Hecke-module for 
many reasons. For example, results of Grunewald and Schwermer imply that, 
for all but a finite set of K, the rational cohomology of SL2^6k) contains cus- 
pidal cohomology, which is important in the theory of automorphic forms. Also, 
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the "Langlands philosophy" predicts a "Shimura-Taniyama-Weil" correspondence be- 
tween Hecke eigenclasses and certain algebraic varieties defined over K. More pre- 
cisely, let r C SL 2 ((?k) be a congruence subgroup and let £ G H 2 (T\$) 3 ; C) be a 
cuspidal Hecke eigenclass. Then one hopes to associate to £ an algebraic variety 
V/K — specifically an elliptic curve or an abelian variety of dimension two — so that 
the zeta function of V is assembled from the eigenvalues of £ in a precise way. Results 
of Cremona 0] and Cremona and Whitley || when K has class number one support 
this. Hence one wishes to compute Hecke eigenvalues for general K. But in general 
the ring &k is not a euclidean domain, and one cannot directly apply the modular 
symbol algorithm as described in [Tj|] (however, see § |1.2| ). 

In this paper we present an analog of the modular symbol algorithm for H d (T; Z), 
where T is a torsion-free arithmetic group associated to certain self-adjoint homoge- 
neous cones, and d is the cohomological dimension of V . This includes finite index 
subgroups of SL 2 (R), where R is 

• Z, 

• the ring of integers in a CM field, or 

• the ring of integers in a totally real field. 

We replace the continued fractions of with a study of the geometry of self-adjoint 
homogeneous cones. Thus our algorithm does not require that R be a euclidean 
domain. 

Here is the organization of this paper. Section contains a review of the classical 
modular symbol algorithm from [18| and presents our algorithm in that case. Section 
|3] contains a review of the reduction theory of self-adjoint homogeneous cones, and 
in Section |] we generalize some results of |22| to this setting (Theorems |2] and |3|). 
These results are valid for cones of any Q-rank > 1. Finally, Section [5] contains a 
description of our algorithm, in Theorem |j. Throughout the paper we comment on 
implementation issues related to the algorithm. 



1.2. Related work. Let T C SL 2 {@k), where 

is the ring of integers in an 
imaginary quadratic extension K = Q(y/—m). For K a non-euclidean ring with class 
number one (m = 19, 43, 67, 163), Whitley developed a "pseudo-euclidean algorithm" 
that allowed implementation of the modular symbol algorithm p3 |. 

Also, I learned upon completion of this work that Jeremy Bygott has independently 
studied the modular symbol algorithm for the non-PID imaginary quadratic case in 
his forthcoming Ph.D. thesis ||. 



1.3. Acknowledgments. The results in this paper depend heavily on results from 
the work of Avner Ash. I thank him for graciously and patiently explaining his work 
to me. I also thank Mark McConnell for a careful reading of an early version of 
this paper and many helpful comments. I also thank the referee for many helpful 
suggestions. 
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Finally, the results in this paper are an extension of some of the results in my Ph.D. 



thesis [14|. I thank heartily my advisor, Robert MacPherson, for the encouragement 



and inspiration he has given me. 

2. A MOTIVATING EXAMPLE 

In this section we illustrate our algorithm for T C SX 2 (Z). As in the introduction let 
9)2 = SL2(R)/ SO2, and let $)% = ^2 U Q U {00} be the usual partial compactification 
of S)2 by adding cusps. We assume that $)% is given the Satake topology, and we 
extend the action of 5 , L 2 (1R) to the cusps. We denote the quotient T\S) 2 by Xr- We 
assume that T is torsion-free, so that Xr is smooth. 



2.1. We begin by paraphrasing aspects of Manin's work |L8[]. By Poincare duality, 
i7 1 (X r ;C) may be identified with ifi(Xr;C), and so we may study the space of 
weight-two modular forms for T by studying the latter. Let q\ and g 2 be cusps 
equivalent modulo V . Then any smooth path 7 from q% to q2 descends to a closed 
path on Xr representing a class in ifi(Xr; Z). Furthermore, this class is independent 
of 7, and in fact depends only on the ordered pair (qi, q 2 ). 

More generally, suppose q\ is not necessarily equivalent to q 2 modulo Y. Then 
integration of one-forms uj G H 1 (X-p] R) along 7 yields a functional J : H 1 (X-p] R) — > 
R, and this allows us to associate to the pair (gi,^) a class in Hi(X r ;R). By the 
theorem of Manin-Drinfeld ([0, p. 61), this class actually lies in if^X^Q). We 
define a modular symbol to be the rational homology class constructed from an ordered 
pair of cusps in this way, and denote this class by [qi, q 2 }. This class agrees with the 
class in the previous paragraph when qi and q2 are equivalent modulo V. 

Proposition 1. |18j The modular symbols satisfy the following: 

1- [9i,?2] = -[92,91]- 

2- [91,92] = [91,93] + [93,92]- 

Furthermore, Hi(Xr] Z) is spanned by modular symbols modulo T. 

2.2. The Hecke operators act on iJ 1 (Xr;C), and hence by duality on Hi(X-p]C). 
On a modular symbol, an operator acts by 



[91,92] 1 — > [«9i,«92], 



where A is a finite set of 2 x 2 integral matrices that depends on the operator. For 
example, let F = T(N), and let p be a prime not dividing N. For the classical operator 
T p we may take 

p Wl OWl l\ / 1 p-l 



" U(7p i 1 j ' 1 p j ' 1 p y ' " ' ' 1 p 
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Here a p G SL 2 (Z) is a fixed matrix satisfying 



P 



-1 







mod N 



20[, Prop. 3.36). Note that deta ^ ±1. 
A finite basis of Hi(Xr\ C) is provided by the set of unimodular symbols. Write a 
cusp q in lowest terms as m/n, where the cusp oo is written formally as 1/0. Then 
the unimodular symbols are the symbols [51,52] satisfying 

detf™ 1 m2 )=±l. 

V n l n 2 J 

The Hecke operators do not preserve unimodularity, and it is necessary for eigen- 
value computations to construct an explicit homology between a non-unimodular 
symbol and a cycle of unimodular symbols. This is done by the modular symbol 
algorithm. Assume that a non-unimodular symbol has the form [0,q], where q is a 
positive rational number. Let [ai, . . . , a J be the simple continued fraction expansion 
of q, i.e. 

1 

q = ai H 



a 2 + 



1 



+ 



OA- 



Let qi be the i th convergent [ai, . . . , a^]. Then by applying (2) from Proposition [1], 
we have 

[0, q\ = [0, 00] + [00, q\\ H h [q k -i, q\- 

Furthermore, the basic properties of simple continued fractions imply that the mod- 
ular symbols on the right are unimodular. Figure [I] illustrates the result for the 
modular symbol [0, 12/5]. 




Figure 1. 12/5 
[2, 5/2] + [5/2, 12/5] 



2 12/5 5/2 

[2, 2, 2] implies [0, 12/5] = [0, 00] + [00, 2] + 



To complete the discussion, we note that ST^Z) acts transitively on the cusps. 
Since the above algorithm is 5 , L 2 (Z)-equivariant, any modular symbol can be written 
as a sum of unimodular symbols. 
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2.3. Now we present our technique for writing a modular symbol as an equivalent 
sum of unimodular symbols. No use will be made of continued fractions; instead, 
we look at the relationship between a geodesic representing a modular symbol and 
a certain tessellation of Sj 2 - in this simple case our algorithm will appear needlessly 
complicated, but it is formulated in a way that will generalize to other settings. It is 
also quite practical for machine computations. 

To begin, we tile Sj 2 with the SL 2 (Z)-translates of the ideal geodesic triangle with 
vertices 0, 1, and 00 (see Figure 0). This tessellation descends to a finite triangulation 
of Xr- The edges of this tessellation are geodesies inducing the unimodular symbols, 
and every unimodular symbol arises in this way. 



00 




1 



Figure 2. Ideal triangles in f) 2 

Given any point x G Sj 2 , let R(x) be the set of vertices of the triangle (or edge) of 
the tessellation meeting x. 

Let [0, q] be a modular symbol as before, and let 7 be the ideal geodesic in fj 2 from 
to q. Because 7 is a geodesic between two rational cusps, one can show that 7 will 
only meet a finite number of triangles in the tessellation. Hence we may choose a 
finite subset x±, . . . , x r G 7 (as in Figure |3|) so that 

LOG 

2. q G R(x r ), and 

3. R{xi) n R{x l+1 ) ^ 0. 

We call such a collection a sufficiently fine partition of 7. From each R(x{) H R(xi + i) 
choose a cusp gj. Then we claim that we have a homology 

(1) M = [0,gi] + [qi,q 2 ] + ••■ + [q r ,q], 

and that each term on the right is a unimodular symbol. 

First, we may see we have a homology either by repeatedly applying (2) of Propo- 
sition H or by continuously deforming 7 into geodesies inducing the classes on the 
right-hand side (see Figure £|) . 

Finally, each nontrivial term on the right of (fl|) corresponds to an edge in the 
tessellation because and q i+ i are both vertices of a triangle containing x i+ \. In 
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Figure 3. A partition of 7 




Figure 4. A deformation of 7 



fact, as shown in [14], if 7 is oriented then one may choose the qi canonically, and our 



algorithm in this case is equivalent to the modular symbol algorithm. 

3. Self-adjoint homogeneous cones 

In this section we present the geometric context of our algorithm. This specifies 
the arithmetic groups to which our algorithm applies and describes the constructions 
replacing S)2 and its tessellation from §0. 

The results in § |3.1| - |3.4] and § |4. 1| are due to A. Ash and originally appeared in JTJ 
and m . Our exposition closely follows the former. 

3.1. Let V be a real vector space defined over Q, and let C C V be an open cone. 
That is, C contains no straight line, and C is closed under homotheties: if x G C 
and A G 1R >0 , then Ax G C. The cone C is called self-adjoint if there exists a scalar 
product ( , ) on C such that 

C = { x G V I (y, x) > for all y G C - {0} } . 

Let G denote the connected component of the identity of the linear automorphism 
group of C, i.e. G = { g G GL{V) \gC = C } u . The cone C is called homogeneous 
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if G acts transitively on C . If K denotes the isotropy group of a fixed point in C, 
then we may identify C with G/K. The fact that C is self-adjoint implies that G is 
reductive and C modulo homotheties is a Riemannian symmetric space. 

We also assume that all these notions are compatible with the Q-structure on V. 
That subgroup of GL(V), G is defined by rational equations, and the scalar 

product ( , ) is defined over Q. This is stronger than saying that G is defined over Q. 
In particular, the group of real points G(M) must be isomorphic to a product of the 
following groups ([^0[, p. 97): 

1. GL n (R) 

2. GL n (£) 

3. GL n (M) 

4. 0(1, n - 1) x R x 

5. The noncompact Lie group with Lie algebra e 6 („ 2 6) © K 

In each case V is a set of hermitian symmetric matrices. In other words, V is the set 
ofnxn matrices over an appropriate IR-algebra with involution r, in which A G V if 
and only if A t = A 7 '. The cone C is then the subset of "positive-definite" matrices in 
an appropriate sense. For details we refer to ( IjlOfl , Ch. V). 

3.2. Let if be a hyperplane in V. We say that if is a supporting hyperplane of C if 
H is rational and H n C = but H n C ^ 0. 

Given a supporting hyperplane H of C, let C" = Int (H fl C). (Here Int(v4) is the 
interior of A in its linear span.) Then C' is called a rational boundary component, 
and is a self-adjoint homogeneous cone of smaller dimension than C. 

Definition 1. The cusps of C are the one-dimensional rational boundary components 
of C. The set of cusps is denoted S(C). 

3.3. Let L C V(Q) be a lattice, i.e. a discrete subgroup of V(Q) such that L ® Q = 
V(Q). Let denote the subgroup of G(Q) preserving L. An arithmetic subgroup 
of G is a discrete subgroup commensurable with Tl for some L. Any torsion- free 
subgroup T C of finite index will act properly discontinuously and freely on 
C. Thus the quotient T\C is an Eilenberg-Mac Lane space for T, and the group 
cohomology H*(T) is H*(T\C). In fact, since homotheties commute with the action 
of T, we may pass to X := M >0 \C, and compute H*(T\X) instead. 

3.4. Let A C V(Q) be a finite set of nonzero points. The closed convex hull o of 
the rays {M-°a; | x E A} is called a rational polyhedral cone. The rays through the 
vertices of the convex hull of A are called the spanning rays of a. We denote the 
set of spanning rays by R(a). The group G(Q) acts naturally on the set of rational 
polyhedral cones, and we denote the action by a dot: o \— > g ■ a. 

Now we would like to partition C into convex subsets using a collection of rational 
polyhedral cones, in a manner compatible with the T^-action. This requires some 
care, as C is open and any a as above is closed. 
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Definition 2. (@, p. 117) Let T C T^ be an arithmetic subgroup of G. A set of 
closed polyhedral cones {c a } is called a T-admissible decomposition of C when the 
following hold: 

1. Each a a is the span of a finite number of rational rays. 

2. For each a, the cone o a C C . 

3. Every face of a a a is some erg in the decomposition. 

4. a a PI erg is a common face of cr Q and erg. 

5. For any cr Q and any 7 <G T, 7<r Q is some erg in the decomposition. 

6. Modulo T, there are only a finite number of cr a 's. 

7. c = U>„nC). 

Note that a T-admissible decomposition descends to a decomposition of X into 
open cells. 

We now describe a technique to construct T-admissible decompositions. The tech- 
nique originates with Voronoi, and was generalized by Ash to all self-adjoint homo- 
geneous cones. Let V be L — {0}. 

Definition 3. The Voronoi polyhedron IT is the closed convex hull of V D H(C). 

Theorem 1. p. 143) The cones over faces ofTl form a T-admissible decompo- 
sition of C . 

We call the cones in this T-admissible decomposition of C the Voronoi cones, the 
decomposition of X associated to IT the Voronoi decomposition, and the open cells in 
X the Voronoi ce//s.Q Here are two examples of this construction. 

Example 1. The original example investigated by Voronoi is the following. Let 
V be the vector space of symmetric n x n matrices, and let C C V be the cone of 
positive-definite matrices. Then G is GL„(IR) , which acts on V by A 1— > gAg 1 , and 
K is SO n . The scalar product is given by (A, B) = Tr(AB). 

Let L be the lattice of integral symmetric matrices. Then T^ = SX n (Z). The 
set of cusps H(C) is obtained as follows. Any nonzero integral (column) n-vector v 
determines a rank one quadratic form by v 1— > w*. The cusps are the rays generated 
by all points in C of this form. Suppose that z G L' H H(C) is a cusp arising in this 
manner from the n-vector t>. Then if y £ C, the scalar product (z,y) is equal to the 
quadratic form y evaluated on v. 

For n = 2 we have that X = Sj2, and the Voronoi decomposition is the tessellation 



described in § [2.3| . For n > 4 there is more than one type of top-dimensional Voronoi 
cone modulo T^, and not all the top-dimensional cones are simplicial. Complete 
information about these decompositions for n = 3 and 4 can be found in [17] and [[19] . 



-'^Our terminology is nonstandard. In |3|, the polyhedron IT is referred to as a "kernel comparable 
to a T-polyhedral cocore" and is a specific example of a more general theory. 
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Example 2. Let K/Q be an imaginary quadratic extension, and let &k be the ring 
of integers of K. Let V be the vector space of 2 x 2 hermitian symmetric matrices 
over C, and let C C V be the cone of positive-definite matrices. Then V is a four- 
dimensional vector space over IR, and X is three-dimensional hyperbolic space Sj^. 
For L we may take the matrices in V with entries in and then = SL 2 (0k). In 
the classical picture of ^ 3 C C x IR-°, the rays generated by the vertices of II become 
the points K U {00}, where K is pictured as a subset of C x {0} and 00 is pictured 
infinitely far above C x {0} along R-°. 

The Voronoi decomposition becomes a tessellation of ^3 into ideal three-polytopes. 
In general these polytopes will not be simplices. For example, if K = Q(y/—1), then 
the unique polytope modulo T L is an octahedron |7[][]I2"|. If &k is euclidean there is 



only one type of top-dimensional polytope in the tessellation, but for general K there 
will be more than one type. 

3.5. Here is the connection between the Voronoi decomposition of X and H*(T; Z). 
Let N be the dimension of X and d the cohomological dimension of V . Let C 
be the set of Voronoi cells of codimension k. The group V acts naturally on C k 
by its action on rational polyhedral cones. We want to construct a r^-equivariant 
"coboundary" map 5 k : Z[C fc ] — > Z[C fc+1 ] so that the resulting chain complex modulo 
T computes if*(T;Z). We will call (C*,S*) a cocell complex, and will say that the 
Voronoi decomposition gives X a cocell structure. 

According to JE], there is a topological space W C X such that the following hold: 

1. W admits the structure of a T-equivariant regular cell complex with top-dimensional 
cells of dimension d. 

2. W is a deformation retract of X, so that the homology of the chain complex 
associated to T\W is the homology of T\X, and hence of T. 

3. This cell structure is dual to the Voronoi decomposition in the following sense: 
every k-ce\\ of W transversely intersects exactly one Voronoi cell of codimen- 
sion k. 

Let Wk denote the set of /c-cells of W. Given r 6 Wk, we denote its dual cell by f G C k . 
Let Z[Wfe] (respectively Z[C fc ]) denote the free abelian group on the elements of Wk 
(resp. C k ). 

We may choose orientations compatibly between the cell and cocell structures in 
the following sense: for each pair (r, f) we may fix orientations so that in the home- 
omorphism 



Int r x f 



the product of the orientations is carried to a fixed orientation of K . This constructs 
a map Z[C7 fc ] -> Hom z (Z[W fc ], Z). 
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Now to construct 5*, we use the boundary map from W*. Given two cells a, r G W*, 
write a < r if o appears in the closure of r. Then dk ■ Z[Wfc] — » Z[Wfc_i] has the form 



^[<7 : r] a, 

IT<T 



where the [a : r] = ±1 keeps track of the relative orientation between a and r. 
(Saying is a regular cell complex makes [a : r] well defined.) We define 5 k : Z[C fe ] — > 
Z[C fe+1 ] by 

(2) f_>£[ T:a ]* 

r<(7 

Proposition 2. W^t/j i/je above coboundary map, H k {T\'L) is naturally isomorphic 
to the k th -cohomology of the quotient modulo F of (Z[C*] , 8*) . 

Proof. We must show that 5 2 = and that 5 is the adjoint of d with respect to the 
pairing between Z[W*] and Z[C*]. The former is purely formal using (0) and the fact 
that f is the map Z[Wjt] — > Z that takes r to 1 and all others to 0. The latter is 
easily verified from the definitions and the choice of orientations. □ 

4. VORONOI REDUCTION 

In this section we address two questions: 

1. How do we construct IT in practice? 

2. Given a point x G C, can we determine a top-dimensional Voronoi cone contain- 
ing it? (Such a cone is unique for generic x, and for any given x there are at 
most a finite number of such cones containing it.) 

In 1 22], Voronoi answers these in the setting of Example |l|, where C is the cone of real 



positive-definite symmetric matrices. In this section we prove that Voronoi's results 
remain true in our more general context. 

4.1. First we describe some geometric properties of II which are proved in [I]. 

Let F be a facet of II, that is, a codimension-one face of II. Then there is a unique 
point y F e C n V(Q) such that 

1. F = { x G II | (x, Up) = 1}, and 

2. for all x G II — F, we have (x, yp) > 1. 

We say that yp defines a supporting hyperplane of II. 

Given a facet F, let Zp be the finite set of points z G L'nH(C) such that (z, yp) = 1. 
Then F is the convex hull of Zp, and as we range over all w G L' H H(C) such that 
u> G" Zjt, the set of numbers (w,yp) is bounded below away from 1. We call yp the 
perfect form associated to F and Zp the set of minimal vectors of yp. In the case of 
Example [l], the are perfect quadratic forms in the classical sense, with minimal 
vectors Z F |22]. 



MODULAR SYMBOLS FOR Q-RANK ONE GROUPS AND VORONOI REDUCTION 



11 



Let a C C be a rational polyhedral cone. Then a satisfies the "property of Siegel" 
with respect to the Voronoi decomposition. Specifically, the intersection a D II is cut 
out from II by a finite number of supporting hyperplanes (0, p. 73). This implies 
that for any a and for any x G C, the orbit T^x meets a in a finite set. 

Given any y G C(Q), let ir(y) : V — ► R denote the linear map x t— > (a;, ?/). We also 
need the following finiteness result. 

Proposition 3. Let y G C(Q). T/ien /or any > ; £/ie sei 

{zGL'nH(C)|0<<2:,y)<//} 

zs finite. 

Proof. Given any A G R, let H\ be the affine hyperplane {x G V | (x, y) = A}. Then 
H = ker(7r) is rational, since y G C(Q) C V(Q). Hence the map 7r(w) takes L onto 
a lattice in R. Since some multiple of w lies in Z/, this lattice is nontrivial. Thus to 
prove the claim, it is enough to show that for any A > 0, the set H\ PI L' n 5(C) is 
finite. 

To see this, consider the set CnH\. Note that H\ meets C, since y' := G H\. 

{y,y} 

Let £ C H\ be any line through y', and let be C \ C . Since C is a cone, £ must 
leave C, and hence £ D 9C7 7^ 0. Let x be a point in £ n The self-adjointness of 
C implies that there is another point z on i D 9c7. (See Figure |^.) 



a; 




z 






j 








3 








Figure 5. 





Now let £ range over all lines in H\ through y', and consider the set of lengths of 
the segments xy' and y'z. This set is bounded above, and hence C fl H\ is a closed 
and bounded subset of V. Since L is a lattice, the set H\ flX'n 5(C), if nonempty, 
is finite. □ 

4.2. Now we describe the construction of IT. Call two facets of II neighbors if they 
meet along a codimension-two face of IT. We show how, given a facet F, one may 
systematically find all neighbors of F. 

Lemma 1. Let F and G be neighboring facets of IT with perfect forms yp and yc- 
Let v G V(Q) be orthogonal to the affine span of the origin and the polytope F nG, 
and such that (x, v) > for all x G F. Then yc = yF + pv for a unique p G M >0 . 
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Proof. First note that the affine span of F D G and the origin is a hyperplane, since 
FC\G is a co dimension-two face of II. Thus v is unique up to a scalar. Let v' = y G —yF- 
Then (x, v') = for all x G F D G, and t> ' 7^ 0. Thus pt> = u' for some nonzero pel, 
which shows that y G = pv + Vf- We must show p is positive. 

Let x G F — (F fl G). Then 1 < (x,y G ) = (x,yp) + p(x,v) = 1 + p(x,v), which 
means that p(x,v) > 0. Since (x,v) > for all x & F, the result follows. 

□ 

Suppose that we are given a facet F with corresponding perfect form yp and min- 
imal vectors Zp. Choose a maximal proper face E C F, and let Zp C Zf be the 
minimal vectors affinely spanning E. Let G be the facet neighboring F along E, and 
write yc = yF + pv, where v is a vector satisfying the conditions of Lemma [I]. Let Z G 
be the set of minimal vectors of y G - Define the function p(x) by 

/ \ 1 - (x,Vf) 

P\X) ■= / v 

(x,v) 

and define S by 

S := {x G L' n E(C) I (x, v) < and yp + p(x)v G C} . 

Note that 5 is nonempty. This follows because whenever z G Z G \ Zp, a computation 
shows p(z) = p and yp + p(z)v = y G , implying z G S. This same computation shows 
conversely that if z G V fl H(G) and p(z) = p, then z is a minimal vector of y G not 
in Z E . 

Lemma 2. As x ranges over S , we have p(x) > p. 

Proof. Assume there is some x' G S with p(x') < p. Then 

= (J,Vf) + p(x',v) 
< (x',y F ) +p(x')(x',v) 
= 1, 

which is a contradiction. □ 

Now we show how to compute y G . Choose any £ G S and consider the point 
Vl '■= Vf + p(£)v. Let T be the set 

T:={xeL'n~(C)\(x,y e )<l}. 
Lemma 3. The set T fl (Z G \ Zp) is finite and nonempty. 

Proof. Since v G V(Q), we have that ye G C(Q). Thus T is finite by Proposition ^. 
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We now show T H (Zg \ Z F ) is nonempty. Let z G [Zq \ Ze) H S, a set we have 
seen is nonempty. Then 

(z,yt) = (z,y F ) + p(£)(z,v) 

< (z, y F ) + p(z, v) [Lemma || 

= (z,Vg) 
= 1. 

Thus z G T, and T H (Z<3 \ Ze) is nonempty. □ 

Theorem 2. Given a facet F of the Voronoi polyhedron, all its neighbors may be 
determined in a finite number of steps. 

Proof. Let y F be the perfect form corresponding to F, and let Z F be the set of minimal 
vectors of y F . Using standard techniques of convex geometry, such as Fourier-Motzkin 
elimination (p4||, p. 37), we may determine all the maximal proper faces of F. Let E 
be such a face, and let Z F C Z F be the minimal vectors affinely spanning E. Using 
Zei we may determine v as in Lemma [1] Let G be the facet of IT neighboring F along 
E, and let yo be the corresponding perfect form. We need to compute ya- 

First we find an i G S by searching over L, and then using I we construct the 
finite set T. The latter can done in a finite number of steps because it is equivalent 
to finding the set of vectors in a lattice on which a positive-definite quadratic form is 
less than a constant (cf. Example [l] and the final paragraph of §3.1| ). By Lemma [3], 
T n {Zq \ Ze) is nonempty, where Zq is the set of minimal vectors of yc- Now let 
Z C (T\ Ze) be the set on which p(x) attains its minimum. By Lemma and the 
paragraph preceding it, Z C Z G \ Z E - Let H be the affine span of Z U Z E - Then yc 
is the unique point satisfying (x, yc) = 1 for all x G H. 

Repeating this procedure for each maximal face of F, we may determine all the 
neighbors of F. □ 

Remark 1. In practice, it may be the case that S only consists of Zq \ Ze-, as a 
computation with SX 2 (Z) shows. 

Hence one may find facets of II provided one can construct an initial facet. In 
the setting of Example |], Voronoi did this by showing that the quadratic form A n , 
defined by 

n 

/ ] X i + ^ ] ( X i ~ ' X j) ' 

i=l l<i<j<n 

is perfect for all n (|[22||, §29) .Q In our more general setting, one cannot write down 
a perfect form that works for every case, even if one restricts to Bianchi groups 
(Example |2|). However, in practice one may do the following. 



2 Voronoi called this form the principal perfect form. 
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First choose a large bounded set U C V containing the origin, and let £ be the 
convex hull of L' fl 5(C) n U. Then £ is a bounded polytope in V. Furthermore, 
since the facets of II are bounded, if U is sufficiently large then most facets of £ will 
be facets of II. 

To check if a facet F C £ is a facet of II, one computes yp and checks whether the 
z £ L' fl 5(C) such that (z, yp) = 1 are vertices of F. 

Remark 2. In the Bianchi case, another possibility is to first construct the retract W 



from § [375| using techniques in |21[ , and then use the duality between them to deduce 
the structure of II. 

4.3. To answer question (2), Voronoi describes a reduction algorithm. This algo- 
rithm is based on the following: 

Proposition 4. Fix an x £ II and a real number \i > 1. JTien t/iere are on/?/ a finite 
number of perfect forms yF satisfying 

(x,Uf) < V- 

Proof. Choose a facet F such that (x,yp) < \i. Since £ C(Q), Proposition |3] 
implies the set £ L' fl 5(C) | (z,yp) < fi} is finite. Hence the polyhedron 

E = {yeU\(y,y F )< f i} 



is a bounded polytope in C. If // £ ir(yp)(L), where ir is defined in § [4.1| , then the 
hyperplanes bounding £ will be rational with respect to the Q-structure on V, and 
thus S will have vertices in V(Q). Hence, replacing \x by a slightly larger number if 
necessary, we may assume £ has vertices in V(Q). 

Therefore the cone generated by the vertices of £ is rational polyhedral, and only 
meets a finite subset of the orbit Tlx. By taking the adjoint action of with 
respect to the scalar product, it follows that only finitely many facets F' that are 
T^-equivalent to F satisfy (x, yp>) < //. Since there are only a finite number of facets 
modulo T L , the result follows. □ 

The following lemma gives a local condition for when a given point of C lies in a 
cone over a facet of n. 

Lemma 4. Let F be a facet of H, and let £f be the set of neighbors of F. Let x £ C. 
Suppose that (x, yp) < (x, ya) for all G £ & . Then x lies in the cone over F . 

Proof. Choose A so that x' := Xx satisfies (x r , yp) = 1 and (x', ya) > 1 for G £ W . 
For e > 0, let £ e be the polyhedron defined by 

£ e = {x | 1 < (x, y F ) < 1 + e and (x, y G ) > 1 for G £ Sf .} 

If e is sufficiently small, then £ e C n. Hence x' £ £ e , and thus x' £ n. Since x' also 
lies in the supporting affine hyperplane {x \ (x, yp) = 1}, it must lie in F. □ 

Now we describe our Voronoi reduction algorithm. 
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Theorem 3. Let x G H, and choose a facet F. Let p, = (x,i/f)- The following 
algorithm determines a cone in the Voronoi decomposition containing x: 

1. For each neighbor G of F, compute (x,yc)- 

2. If there exists a neighbor G with (x,yc) < p, replace F with G, p with (x,yc), 
and return to step one. 

3. Otherwise, terminate the procedure: x lies in the cone generated by F . 

Proof. By Lemma U if the algorithm terminates then we have determined a cone 
containing x. We now prove that the algorithm terminates. Suppose that a neighbor 
G of F has (x, yc) < P- (Note that this quantity must be positive since C is self- 
adjoint and x, y G C.) Then we return to step one, and we have decreased the scalar 
product. Since by Proposition [| the set of facets satisfying (x,ypi) < p is finite, the 
algorithm must terminate. □ 

Remark 3. The data needed to implement this algorithm is the same as that needed 
for the structure of II modulo T^, along with some additional information. In partic- 
ular, one must determine: 

1. A finite set of representatives of the facets of II modulo r^. 

2. For each F G and each neighbor G of F, a element 7 G Tl such that 7G G & . 
For an example of the implementation of this algorithm for SX2(Z), we refer to [[15 



As far as we know, the computational complexity of this algorithm is unknown. 
However, in our experience it performs very well for SX 3 (Z) and SX 4 (Z). 

5. The modular symbol algorithm 

In this section we define modular symbols and describe our Hecke algorithm. Our 
definition of a modular symbol is closely related to the definitions appearing in || 
and 0). As before, let N be the dimension of X and let d be the cohomological 
dimension of V. We assume from now on that G has Q-rank one, so that N = d + 1. 

5.1. Let p(II) be the set of rays in C generated by the vertices of the Voronoi 
polyhedron II. Given an ordered pair (u,v) G p(II) x p(II), we want to construct a 
class [u,v] G H d (T;Z). 

To this end we recall that X may be extended to a bordification X such that the 
quotient T\X, the Borel-Serre compactification, is a compact manifold with corners 
with interior T\X 0. Let 7r: X — > T\X be the canonical projection. Given u,v G 
p(II), we determine a path in T\X as follows. First u and v determine a closed cone 
a C V, and we let a be the closure of a fl V mod homotheties in X. Then ir(a) is a 
path with endpoints lying in d(T\X). Choosing an ordering (u, v) fixes an orientation 
of 7r((j) and thus determines a class in Hi(T\X, <9(r\X); Z). By Lefschetz duality, 

tfa(r\Z, d(T\X); Z) = H d (T\X; Z), 

and since T\X is homotopy equivalent to T\X, we have actually determined a class 
mH d {T;Z). 
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Definition 4. A modular symbol is a class in H d (T; Z) constructed as above from an 
ordered pair (u,v) G p(jT) x p(IT). The class is denoted [u,v]. 

Note that this definition is almost the same as that in for T C SX 2 (Z), because 
using (u, v) to generate a cone is the same as choosing a specific path as in § [2.1| . 
However, as in the earlier definition, the class [u, v] is independent of this path. 

We have the following analogue of Proposition [[]: 

Proposition 5. Let u,v e p(n). T/ie modular symbols satisfy the following: 

1. = -[v,u). 

2. If w E p(II) ; t/ien [it, = [it, w] + [w,v]. 

3. The modular symbols span H d (T; 7*). 

Proof. Only (2) and (3) require proof. To prove (2), let a be the cone generated by u 
and v, and choose a ray a; C a distinct from u and i> . Let : [0, 1] — > C be a continuous 
family of rays such that 0(0) = x and 0(1) = w. Let c\{t) (respectively cr 2 (i)) be 
the cone generated by u and 0(t) (resp. 0(t) and v ). Then provides a continuous 
deformation of a into cxi(l) U 02(1) that induces the homology [u, v] = [u, w] + [w, v}. 

Now to prove (3), note that the results of § [3.5| imply that any class in H d {T]'L) 
may be written as a cocycle for 5 using the Voronoi cones of co dimension- d modulo 
T. But the Voronoi cones of codimension d are cones generated by pairs of vertices 
of n, and so any such cycle is in the span of the modular symbols. □ 

Remark 4. The third statement is also a consequence of the results in 0. 

5.2. Now we describe our replacement for the notion of a unimodular symbol. Recall 
that R(cr) is the set of spanning rays for a cone a (§ |3.4| ). 

Definition 5. Let a C C be a rational polyhedral cone. Then a is Voronoi-reduced 
if there is a top-dimensional Voronoi cone a a such that 

R{a) C R{a a ). 

Note that a Voronoi-reduced cone is spanned by cusps. A Voronoi-reduced cone 
need not be a Voronoi cone, since the top-dimensional Voronoi cones are not simplicial 
in general (cf . Examples |I] and |2|) . However, since every Voronoi cone is generated 
by a finite set of cusps, and because of the finiteness properties of any T-admissible 
decomposition, we have 

Proposition 6. Modulo V , there are only finitely many Voronoi-reduced cones. 

We say that a modular symbol is Voronoi-reduced if it is induced by a Voronoi- 
reduced cone. Propositions |5| and pimply that the Voronoi-reduced modular symbols 
provide a finite spanning set for H d (T] Z). 

Remark 5. Although the Voronoi-reduced modular symbols provide a finite spanning 
set for H d (T] Z), they are not a basis of H d (T; Z). In fact, they are not even a basis of 
Z[C*] modulo T, because they are not necessarily supported on codimension-(i Voronoi 
cones. However, in practice this does not affect their usefulness (cf. Remark ^). 
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5.3. Let u,v G p(IT), and let a be the rational polyhedral cone generated by u and 
v. We are now ready to describe and prove our algorithm. 

Given x G C, let a(x) be the unique Voronoi cone containing x, and by abuse of 
notation let R(x) be R(a(x)). 

Let xi, . . . ,x r be points in a fl C such that the rays M >0 xi, . . . , M >0 x r are distinct. 
These points subdivide a into a collection of cones, namely those generated by the 
pairs (u, xi), . . . , (x r , v). 

Definition 6. The decomposition of cr by the Xi as above is called a sufficiently fine 
partition if 

1. ii€ R(xi), 

2. v G R(x r ), and 

3. for i = 1, . . . , r — 1, we have i?(xj) fl i?(xj + i) 7^ 0. 
Lemma 5. Sufficiently fine partitions 0/7 exist. 

Proof. Since <r is rational polyhedral, the Siegel property implies that a fl IT is cut out 
by finitely many supporting hyperplanes. Hence a fl C meets only finitely many top- 
dimensional cones, and the intersection of these cones with a subdivides the latter 
into finitely many 2-cones {o"i}. We may take Xj to be any nonzero point in the 
interior of ov Conditions (1) and (2) of Definition |6| are trivially satisfied. Also, 
R(xi) fl R(xi + i) 7^ because o~i C Vi and a i+ i C Vi+i, where Vi and Vi+i are Voronoi 
cones that have a face in common. 

□ 

Remark 6. For computational purposes, one may construct sufficiently fine partitions 
of 7 as follows. Let it, v G p(n) and let a be the cone generated by u and v. Choose 
points u G u and v E v. Let x be the midpoint of the segment between u and v, and 
let x be the cone generated by x. Now apply Theorem [I] to check whether u G 
and v G R(x). If these conditions are not satisfied, bisect the segments between u, x 
and x, v, and check conditions (1), (2), and (3). Eventually, by the Siegel property, 
after a finite number of iterations one will have constructed a sufficiently fine partition 
of a. 

Now we present our algorithm. 

Theorem 4. Given a modular symbol [u,v], the following constructs a chain of 
Voronoi-reduced modular symbols homologous to [u, v] : 

1. Choose a set of points {x^ inducing a sufficiently fine partition of the cone 
generated by u and v. 

2. For i — 1, . . . ,r — 1, choose a ray qi G R{xi) fl R(x i+ i). 

Then [u,v] = [u,qi] + [qi,q 2 ] H h [q r ,v\. 
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Proof. First note that each modular symbol on the right hand side is Voronoi-reduced, 
since qi and gj+i are both rays from R(xi + i). We must show there is a homology 
between the right side and the left. Notice that 

[u,v] = [u,qi] + [qi,v] 

by Proposition ||. Repeatedly applying this proposition, we see that 

fa,v] = [qi,qi+i] + [qi+i,v], 

which completes the proof. □ 

Remark 7. For computational purposes, to determine the action of a Hecke operator 
we must write any modular symbol in terms of a basis of H d (T; Z), and by Remark |5| 
the technique in Theorem ^ is not sufficient to do this. However, in practice we 
may precompute explicit homologies between Voronoi-reduced modular symbols and 
modular symbols supported on Voronoi cones, as follows. 

Let & be a set of representatives of the facets of II modulo T^. For each F G J^, let 
u, v be any two vertices of F. Then [u, v] is a Voronoi-reduced modular symbol. To 
write [u, v] in terms of the basis of codimension-rf Voronoi cones, choose any sequence 
of vertices u — u , ui, . . . , = v of F such that Ui and u i+ i are joined by an edge of 
F. Then 

[u,v] — [u, Ui) H h [ttfc_i,v] 

is the desired homology. Now repeat for all u, v and all F G & . 
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